Integrate[Integrate[Integrate[-q*sin[t]*cos[p]*(R+q*cos[t])*cos[p]*sin[t],{t,0, 2Pi}] ,{p,0, 2Pi}] ,{q,0, r}